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Abstract 

Spirals have been studied from a dynamical system perspective starting with Barkley's 
seminal papers linking a wide class of spiral wave dynamics to the Euclidean symmetry of 
the excitable media in which they are observed. However, in order to explain certain non- 
Euclidean phenomena, such as anchoring and epicyclic drifting, LeBlanc and Wulff introduced 
a single translational symmetry-breaking perturbation to the center bundle equation and 
showed that rotating waves may be attracted to a non-trivial solution manifold and travel 
epicyclically around the perturbation center. 

In this paper, we continue the (model-independent) investigation of the effects of inho- 
mogeneities on spiral wave dynamics by studying epicyclic drifting in the presence of: a) 
n simultaneous translational symmetry-breaking terms, with n > 1, and b) a combination 
of a single rotational symmetry-breaking term and a single translational symmetry-breaking 
term. These types of forced Euclidean symmetry-breaking provide a much more realistic 
model of certain excitable media such as cardiac tissue. However, the main theoretical tool 
used by LeBlanc and Wulff can only be applied to their particular perturbation: we show 
how an averaging theorem of Hale can be modified to analyze our two more general scenarios 
and state the conditions under which epicyclic drifting takes place in the general case. In 
the process, we recover LeBlanc and Wulff 's specific result. Finally, we illustrate our results 
with the help of a simple numerical simulation of a modified bidomain model. 
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1 Introduction 

Spiral are found in numerous excitable media [1,2,4,9,13,16,21,22,25,27,35,37,39,40] 
and they give rise to beautiful imagery. While this in itself might yield enough interest 
to study them, there is also (at least) one serious reason to do so: spiral waves have 
been linked to cardiac arrhythmias (to disruptions of the heart's normal electrical 
cycle) [9,17, 36, 37] . Furthermore, 

most arrhythmias are harmless but if they are "re-entrant in nature and [...] 
occur because of the spatial distribution of cardiac tissue" they can seriously 
hamper the pumping mechanism of the heart and lead to death [17, p. 401]. 

As a result, a fuller understanding of spiral wave dynamics in these media becomes 
imperative. 

The equivariant dynamical system approach 

In recent years, one of the most rewarding approach to the study of spiral waves is 
based on Barkley's initial observation that the observed transition from rotating to 
modulated rotating wave can be explained via a Hopf bifurcation together with the 
underlying Euclidean symmetries of the governing reaction-diffusion equations [1,2] 
(i.e.: the semi-flow generated by the dynamical system commutes with the 

u(t, x) i — > u(t, xi cos 9 — x 2 sin 9 + Pi, X\ sin 9 + x 2 cos 9 + p 2 ), (1.1) 

where (9,p 1 ,p 2 ) e§'xR 2 ~ SE(2) and x G M 2 [10,38]). This lead Barkley to formulate 
a simpler ad hoc 5-dimensional ODE system with Euclidean symmetry replicating the 
above transition [3]. 

Sandstede, Scheel and Wulff then proved a general center manifold reduction theo- 
rem (CMRT) for relative equilibria and relative periodic solutions in spatially extended 
infinite-dimensional Euclidean-equivariant dynamical systems, providing a mathemati- 
cal justification of Barkley's insight [11,30-33]. However, this center manifold reduction 
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theorem requires that the spiral wave satisfy certain spectral gap conditions, which of- 
ten fail [34]. Sandstede and Scheel have developed a comprehensive theory of spiral 
instabilities using techniques of spatial dynamics [28,29] to deal with such a situation. 

Other methods are also used to reduce the dynamics to finite-dimensional systems 
(such as the kinematic model using the curvature of the wave as a driving mecha- 
nism [24]), but the equivariant dynamical system approach has the advantage that it 
can often provide universal, model-independent explanations and predictions regard- 
ing the dynamics and bifurcations of spiral waves. For example, the fore-mentioned 
'Hopf bifurcation' from rigid rotation to quasi-periodic meandering has been observed 
in numerically [4] and experimentally [21]. Another example is provided by the anchor- 
ing/repelling of spiral waves on/from a site of inhomogeneity, which has been observed 
in numerical integrations of an Oregonator system [25], in photo-sensitive chemical 
reactions [40] and in cardiac tissue [9]: using a model- independent approach based 
on forced symmetry-breaking, LeBlanc and Wulff showed that anchoring/repelling of 
rotating waves is a generic property of systems in which the translation symmetry is 
broken by a small perturbation [19] J* In the same vein, certain dynamics of spiral 
waves in anisotropic media, such as phase-locking and linear drifting of meandering 
spiral, have been shown to be generic consequences of rotational symmetry-breaking 
[18,26,27]. 



The basic viewpoint 

Consider a piece of cardiac tissue on which numerous (roughly) circular ablation have 
been performed, perhaps in order to treat a patient who is suffering from atrial fibrilla- 
tion [14,23]. These surgical procedures affects both the geometry and the excitability 
of the tissue. Under certain modeling assumptions, any system used to model the 
electrical activity of the tissue needs to incorporate translational symmetry-breaking 
(TSB) components to model the effects of the circular ablations, and a rotational 
symmetry-breaking (RSB) component to model the effects of anisotropy. Let us model 
the electrical properties of such a perturbed piece of anisotropic cardiac tissue using 

^dn this paper, we use the terms 'generic' and 'typical' interchangeably: the set of coefficient values for which 
the anchoring/ repelling property fails to hold has measure zero in the complete coefficient space. 
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a modified version of the bidomain equations of cardiology, under the modeling as- 
sumption that the circular ablation (inhomogeneous) zones consist of a finite number 
of independent "sources" which are localized near distinct sites (j, . . . , ( n in the plane 
(see [7] for a similar hypothesis). The model then has the form 

I 3 n 

«t = ~(U - y - V) + V 2 U + 1 + a ^_ e N **i*i • " OH 2 , M) 

- (12) 
v t = q(u + 6-jv) +22fij9j(\\x-Cj\\ 2 

i=i 

V 2 ^ + e0(a,e)* a!aa!a = eh(a,e)u X2X2 , 

where u is a transmembrane potential, u controls the recovery of the action potential, 
is an auxiliary potential (without obvious physical interpretation), x\ is the preferred 
direction in physical space in which tissue fibers align, e is a measure of that prefer- 
ence, g and h are appropriate model functions, a,s, S and 7 are model parameters, 
/i = (/ii, . . . ,/i n ) G M" is a small parameter and g^ ,v are smooth functions, uniformly 
bounded in their variables [5,8,20,26]. 

If the tissue has equal anisotropy ratios (i.e. e — 0) and the inhomogeneities have 
no effect on spiral wave dynamics (i.e. \i — 0), (11. 2p decouples into the FitzHugh- 
Nagumo equations for u and v, and Poisson's equation for \l/ [26]. 

Let §E(2) denote the group of all planar translations and rotations, and fix an in- 
teger 1 < f and £ G M 2 . The subgroups Zj*+M 2 (the notation will be explained later) 
and SO(2)^ of §E(2) consist of all cartesian pairs of translations and rotations about 
the origin by an integer multiple of 2ir / ' f radians, and of all rotations about the point 
C, respectively. Let V = C+Z^. or T = SO(2)^. Then, T < §E(2) and we will say that 
the semi-flow ^t,e,ti is T— equivariant if it commutes with the restriction of (II. lft to T. 

In the equivariant dynamical system approach, the particular form of the functions 
g™' v is unimportant; the analysis is driven by the fact that (\1.2\i can sustain spiral 
wave propagation [5,8,20,26] and by the equivariance properties of the semi-flow $t i£iA1 
generated by (11.21) . namely: if we neglect boundary effects, the semi-flow 

(El) is SE(2)— equivariant when (e,fi) = 0; 
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(E2) is Z 2 +R 2 — equivariant when e 7^ is small and \x = 0; 

(E3) preserves rotations around Q (but generically not translations) when e = and 
fij = for all j ^ i, and 

(E4) is (generically) trivially equivariant when (e, fi) is a generic small parameter 
vector. 

This is but a special case of a more general family of semi-flows for which (E2) is 
replaced by the following property: the semi-flow 

(E2') is Zj.-i-M 2 — equivariant when e 7^ is small and fi = for some integer j* > 1. 

In [5-7], we used the dynamical system approach to study spiral anchoring in media 
satisfying (El), (E2'), (E3) and (E4): the predictive power of the method was used 
to show that in the case n > 1, spiral anchoring typically takes place away from the 
inhomogeneities. At the time, such a statement defied experimental wisdom. 

Epicyclic drifting 

At this stage, nothing has been said about the nominal topic of this paper: epicyclic 
drifting. The various spiral motions observed in experiments and simulations have 
been classified according to their tip path, an arbitrary point on the wave front that is 
followed in time [9,21]: for instance, the tip path of a (rigidly) rotating wave is a perfect 
circle. Barkley [4] and Wulff [38] have shown that the appearance of an epicyclic tip 
path can be linked to a 'symmetric Hopf bifurcation:' when that happens, every spiral 
wave in the excitable medium is epicyclic. 

However, other epicyclic behaviour cannot be explained by this mechanism. When 
the sizes of the physical domain and of the spiral core are comparable, the latter 
is sometimes attracted to the boundary of the domain and rotates around it in a 
meandering fashion. This has been observed in experiments and numerical simulations 
in a light-sensitive BZ reaction [39,40]. 

Yet another instance of epicyclic motion is shown figure 1: in a bounded region, 
all solutions are attracted/repelled to/by an epicyclic solution manifold. This type of 
spiral wave motion is what we refer to as epicyclic drifting. A la Poincare-Bendixson, 
if a system has a stable epicyclic manifolds (stable in the sense of Lyapunov) it will 
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Figure 1: Epicyclic motion on the stable epicyclic manifold 8. The arrows indicate the direction 
of the flow, while S corresponds to a repelling (perturbed) rotating wave solution pinned at the 
inhomogeneity indicated by the black dot (see [7] for details on spiral anchoring). 



also have a repelling rotating wave (see figure 1), and vice-versa. As such, these man- 
ifolds cannot be observed in fully Euclidean media. What then, can forced Euclidean 
symmetry-breaking (FESB) tell us about epicyclic drifting in systems with the equiv- 
ariance properties of (El), (E2'), (E3) and (E4)? The only work in this vein has 
been performed by LeBlanc and Wulff in [19], in the case n = 1 and without rotational 
symmetry-breaking: unfortunately, the main tool in their analysis cannot be used in 
the general case. 

Article Overview 

The main object of analysis in the present paper is a finite-dimensional system of 
ODE that share the equivariance properties of (El), (E2'), (E3) and (E4) when 
n > 1: it is derived in section 12.11 Then, in section I2.3[ we present a preliminary 
result about averaging which will subsequently be used to prove our main results: to 
wit, when e = or f = 1 and certain conditions are satisfied, there is a (minimal) 
parameter wedge region in which an epicyclic manifold persists. In the case j* > 1, the 
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Figure 2: On the left, the epicyclic parameter wedge regions for the case n = 2 without anisotropy. 
On the right, the epicyclic deleted neighbourhood in parameter space for the case n = 1 with 
anisotropy characterized by j* > 1. If the semi- flow has an epicyclic solution manifold for a 
particular set of parameter values, then the semi-flow has an epicyclic solution manifold (of the 
same stability type) for all parameter values in the adjacent region. Note that these manifolds 
continuously deform along any path contained entirely in the parameter region. The local analysis 
does not provide a clear picture of the behaviour as a path leaves a parameter region. In [5], for 
instance, we give an example where the epicyclic solution manifold disappears as a result of a 
saddle-node bifurcation of rotating waves. The notation will be explained later in the paper. 



epicyclic manifold persists in a deleted neighbourhood of the origin. The parameter 
wedges are illustrated in figure 2, for the case n = 2. The conditions needed depend 
on the kind of forced symmetry-breaking under consideration: in section 3, we study 
the general semi-flow under n TSB terms (i.e. n > 1, e = 0); in section 4, we study 
the combination of a single RSB and a single TSB term (i.e. n — 1, e ^ 0). We 
then combine these results in section 5 to obtain the epicyclic drifting theorems under 
general FESB. Next, we perform a simple numerical experiment on 11.21) with n = 1 
showing the predicted epicyclic motion: to the best of our knowledge, the figure in 
section 6 is the first observed instance of epicyclic drifting in a numerical simulation 
of excitable media. Finally, we give the proofs of two technical results in appendix A. 
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2 Preliminaries 



We start with a derivation of the appropriate center bundle equations describing the 
essential dynamics of spiral waves near a rotating wave under full Euclidean symmetry- 
breaking (FESB). More details on these manipulations can be found in [5,7,30,31,38]. 
Then, in order to lighten the text, we introduce some necessary definitions. Finally, 
we state an averaging theorem which will be used in later sections of this work. 

2.1 The Center Bundle Equations 

In order to facilitate the subsequent analysis, we make the same simplifying assump- 
tions and adopt the same notation as in [6,7]. 

In particular, let X be a Banach space, W C I x K" a neighborhood of the origin, 
and &t,e,fi be a smoothly parameterized family (parameterized by (e, /i) £ U) of smooth 
local semi-flows on X, and let 

a : SE(2) — ► GL(X) (2.1) 

be a faithful and isometric representation of SE(2) in the space of bounded, invertible 
linear operators on X. For example, if X is a space of functions with planar domain, 
a typical SE(2) action (such as ( 11. ip in the preceding section) is given by 

(a( 1 )u)(x)=u( 1 - 1 (x)), 7 e§E(2). 

In this paper, we concern ourselves with the study of epicycle drifting in the case where 
the two following hypotheses are satisfied. The first one is simply a re-telling of the 
equivariance properties (El), (E2'), (E3) and (E4), while the second postulates the 
existence of a rotating wave in the unperturbed SE(2)— equivariant semi- flow. 

Hypothesis 1 There exists 1 < f £ N, distinct points Ci> • • • > d in M. 2 such that if Cj 
denotes the j th vector of the canonical basis in M. n , then Wu £ X, e ^ 0, q^O, t > 0, 

$t,o,a^(a(7)«) = a(7)$t.o,a^(«) 7^ SO(2) Cj , and 

$t,o,o(a(7)w) = a(7)$ t , 0>0 (u), V 7 £ §E(2). 
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Hypothesis 2 There exists u* G X (with trivial isotropy subgroup) and Q* in the Lie 
algebra o/SE(2) such that e n * 1 is a rotation and $ f:0 (w*) = a(e nH )u* for all t. More- 
over, the set {fj,eC\ > 1 } is a spectral set for the linearization a(e~ n *)D§i } o(u*) 
with projection such that the generalized eigenspace range(P*) is three dimensional. 

As discussed previously, such semi- flows can arise from the family of perturbed reaction- 
diffusion systems (1.4) from [7] if e — 0, as well as from the modified bidomain model 
(11.21) given in the introduction if f = 2, for instance. 

It has been shown in [5-7] that, for small parameter vectors (e, /x) G R x M n , the 
essential dynamics of the semi-flow &t,e,n near a (hyperbolic) rotating wave is (locally) 
equivalent to the semi-flow of the following ordinary differential equations on the bundle 
C x S 1 : 



if 

p = e 



+ 0G(t, + V*i((P - (P - A) 



(2.2) 



where v G C, (/3, A) G R x M n , ^,...,( b gC are all distinct and the functions G, Hj 
are smooth, periodic in t and uniformly bounded in p, and G is 2n/j*— periodic in t. 
The specific form of the perturbations is a consequence of forced Euclidean symmetry- 
breaking from the SE(2)— equivariance of (12. 2ft under the following §E(2)— action on 
the bundle CxS 1 : 

(x,6)-(j>,<p) = (e i0 p + x t <p + e), (2.3) 

for all (p,<p) G C x S 1 and (x,9) G SE(2) = C x S 1 , where SE(2) = C+S 1 with 
multiplication (pi, fi)-(p2, ^2) = {e % ' pl P2+Pi, ^1+^2)- The non-standard multiplication 
is made explicit by using the semi-direct product notation +. 

Let f > 1 be an integer and (6 C. In §E(2) = C+S 1 , the subgroup of rotations 
around £ is given by 

§J = {(£,o)-(o, £)•(-£, O^flGS 1 }, 

while the subgroup containing all translations and rotations by angle ^jr, k G Z, is 

C+Zj» = ( ( x, — J :fcGZ:i£C 
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Then C+Z f ~ Z f +R 2 and §^ ~ SO(2) c : under the action described by Q, the 
center bundle equation (12.21) is 

(CI) C+S 1 — equivariant when A) = 0; 

(C2') C+Zj* —equivariant when (3 ^ is small and A = 0; 

(C3) S| £ —equivariant when /3 = and Xj = for all j ^ £, and 

(C4) (generically) trivially equivariant when {(3, A) is a generic small parameter vector. 

Clearly, (12.21) shares the equivariance properties of Hypothesis 1. As such, G 'models' 
the RSB perturbation while the various Hj 'model' the various TSB terms. 

It might seem strange that the parameters (e, fi) are replaced by A) in (12. 2p . 
just as the (j G R 2 are replaced by £ 3 - G C, but since the center manifold reduction 
theorems of [30-33] do not provide an explicit relation between the coefficients of the 
original system of partial differential equations and the reduced ordinary differential 
system of center bundle equations, one cannot conclude that the parameters are the 
same in both systems. 

2.2 Definitions 

An integral manifold is stable if it has a neighbourhood in which all originating positive- 
time solutions approach the manifold exponentially; it is hyperbolic if the linearization 
of the flow on this manifold admits no critical eigenvalues. Let > 0, A > 0, V C R p , 
S = R x V x [0, oq], f : S -> W, g : V -> and h : R x V -> M 9 . We say that / is 
Lipschitz in Hale's sense, which we denote by / G Lip(x; E, i](a, V)), if / is continuous 
in all of its arguments and is Lipschitz in x for (i, x, a) G S with continuous Lipschitz 
constant. 

Next we say that g is bounded by A over V, which we denote by g G B(A;V), 
if ||<7(2;) || < A for all x G V. Finally, we denote the fact that h is T— periodic in 
(p G R by h G ^P^. When the sets S and V are understood from the context, they are 
omitted. Finally, by abuse of notation, we shall often denote O (\x\\ + ■ ■ ■ + \x m \) by 
0(xi, . . .,x m ). 
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2.3 A Generalized Averaging Theorem 

Averaging methods are used to determine whether a particular system has an non- 
trivial invariant integral manifold by studying an averaged system. The main theorem 
is a modified version of one of Hale's averaging theorems (see [5] for details); it can 
easily be extended to the case where v is a parameter vector in W 1 . 

Theorem 2.1 (modified from [15], theorem 6.1, pp. 526 — 527) Let cr > 0. Con- 
sider the system of equations 

(2.4) 

-0 = d(e, is) + 6(t, -0, x, e, v), 

where ip G M, x G [— cr , oo), r y eM ^ depends continuously on (e, v), and d is defined 
over S = [— e , e ] x [— u , v ], with d(0, 0) = 1. For a > 0, let 

E„ = RxRx [-a, a] x S and S = R x R x {0} x S . 

Suppose 9,Ae^n^ and that 
(i) G and A are real-valued over S CT0 ; 
(ii) 6, A G B(E(e, u); S ) where S(e, z/) = 0(e, z/)/ 

(m) /or a// < a < a , G Lip(?/>, x; £ CT , ^(e, u, a)) and A G Lip(ip, x; E CT , ?7(e, z/, cr)), 
wzi/i #(e, z/, cr) = 0(e, z/, a) and r/(e, z/, cr) = 0(e, u, a). 

Then, there exists (ei,z/i) G (0,eo] x (0, z/o] swc/i that for all 

(e,u) G Si = [-ei,e!] x [-1/1,1/1] 
mf/i e ^ 0, (|£.^p /ias a hyperbolic integral manifold T ev which can be represented as an 



invariant torus x = T eiU (t,?p), where 



T e>u G B(D{e, u)) n Lip(0, fi(e, 1/)) n n qr 



■ 



with D(e,u),Q(e,u) —> uniformly as (e, z/) — > 0. Furthermore, the stability of T e u is 
exactly determined by the sign 0/670,0- 
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3 Epicyclic Drifting For n Simultaneous TSB Terms 



When j3 = 0, (12.21) gives the dynamics near a hyperbolic rotating wave for a parame- 
terized family of semi-flows <&t,o,\ satisfying the forced-symmetry breaking conditions 
in hypothesis [TJ We start with a brief review of epicyclic drifting in the case n = 1, 
which was studied in detail in [19], and then present our new results in the general 
case n > 1. 



3.1 The Case n = 1 

Without loss of generality, we may assume £i = 0. In this case, the center bundle 
equation (12. 2p reduce to 

p = e u [v + X 1 H 1 (pe- u ,pe u ,(3)}, (3.1) 
where v G C x and Ai 6 M. is small. Set H(w,w, Ai) = Hi(w — iv,w + iv, Ai). 
Theorem 3.1 ([19], re — written to fit the current symbolism) Let 



I(p) = Re 



/ e~ lt H (pe~ l \ pe u , 0) dt 
Jo 



If po > is a hyperbolic solution of I(p) = 0, then for all Ai 7^ small enough, the 
center bundle equation {\3.1\i has an integral ( solution ) manifold £\ around the origin, 
whose stability is exactly determined by the sign of \\I\po). 

These solutions represent quasi-periodic motion around the origin in the p— plane and 
are observable as epicycle-like motion along a circular boundary in the physical space, 
with angular frequency 1 + O(Ai). Note that the hypotheses of theorem 13.11 are not 
generic: in a random system, I(p) may very well not have a positive hyperbolic root. 

The presence of a repelling integral manifold could explain the fact that spirals are 
sometimes observed to be repulsed by an inhomogeneity if the spiral tip is located 
beyond a certain distance from the perturbation center [25]. 
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3.2 The Case n > 1 

However, the main averaging tool used in [19] to obtained theorem (13. 1 p cannot be 
used to analyze the situation in the case n > 1; furthermore, this difficulty yields a 
interesting twist, as we shall see in this section. 

By re-labeling the indices in (I2.2p if necessary, we can temporarily shift our point 
of view so that £i plays the central role in the following analysis. Set Hj = £j — £i for 
j = 1, . . . , n. Then, under the co-rotating frame of reference z = p — £i + ie lt v , (12.21) 
becomes 

n 

z = e u J2 X j H j{( z ~ -j) e ~ U - iv , ( z ~ Zj)e u + iv, A) . (3.2) 

3=1 

When Ai 7^ and A2 = • • • = A n = 0, we find ourselves in the situation described in 
the previous subsection. Now, set e = Ai, V\ — 1 and Xj = vje for j = 2, . . . ,n, and 
v = (z/2, . . . , u n ) G IR™" 1 . Then (13.21) can be viewed as a perturbation of the corre- 
sponding equation in the case n — 1. Note that S x = and A = (1, v)e. 

Equation (13.21) rewrites as 

n 

z = ee u "j H j{( z - E j) e ~ U - iu » ( z ~ ™> (!) u ) e ) ■ (3-3) 

3=1 

Let Hj(w, w, e, v) = Hj (w — iv,w + iv, (1, v)e) for j = 1, . . . n. Then (13.31) becomes 

z = ee u K(ze- u , ze" , t, e, v) (3.4) 

n 

where K{w, w, t, e, v) = VjHj{w — H^e -4 ', w — S^e**, e, v) is 27r— periodic in t. Con- 

3=1 

sider the near-identity change of variables 

z = w + €k(w, w, t, e, v) (3.5) 
where k G ^ is differentiable in all of its variables. Then 

z = w + e — + — iw + — ty I . 
\ at ow ow 
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Introducing the equivalent complex conjugate equation, this last system becomes 



Kit i) tv 



W 
W 



Z )-e Kt 

Z \K t 



(3.6) 



where k w , n t , k w , k w , n t are used to denote the partial derivatives of k and k. Set 



X = I 2 + e 



Combining (13.61) with (13.41) yields 



w 



w 



A ( e lt K({w + 6K)e lt ,(w + €K,)e lt ,t,e,v) — n t 
e~ lt K ((w + en)e~ lt , (w + eK)e lt , t, e, v) — H t 



(3.7) 



By Taylor's theorem, there are appropriate continuous bounded functions A±, A 2 and 
A 3 G y$ 2 t K satisfying 

e il K ((w + eK)e"*', (w + e«)e", t, e, v) = e u K (we"", fe", t, 0, v) + eAi(w, W, t, e, i/) 

ftt(u>, w, t, e, ^) = Kt(u7, uJ, t, 0, ^ ) + ev4 2 (u>, w, t, e, ^ ) 



and 



where 



T 



I 1 eK TU 
6/^2, 1 



e A 3 (w,w,t,e,u), 



— K w(w, w, t, 0, v) and «^ = k w (w, w, t, 0, v) 



With these, (13 .7p re- writes (upon dropping the equivalent complex conjugate equation) 
as 

w = e (e^K^we^, HJe**, t, 0, v) - K t (w, w, t, 0, u)) + e 2 H(w, w, t, e, u), (3.8) 
where 7i G is bounded and continuous in all its variables. Denote the average 

it wn^it 



value of e lt K(we lt , we l \ t, 0, v) by 



2tt 



M 1 (w,w,z/) = — / e a K(we- it ,we tt ,t,0,v)dt 
2vr y 



(3.9) 
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Then 

e^iue - * We ft ,t,0,i/) = M 1 (w,w, v) + F{w,w,t,u), 
where F G tyf* is uniformly continuous and 



/ F(w, w, t, v) dt = 0. 



(3.10) 



Jo 



Let k be an antiderivative of F with respect to t. Then k e ^Jf^ by (13.101) and 



It is easy to see that M 1 (w,t<;,0) is S 1 — equivariant (see appendix for details); as 
such, there is a continuous function L\ : R — >• C such that M 1 (w,mJ, 0) = wLi(ww) 
[12, p. 360]. 

By Taylor's theorem, there are appropriate continuous bounded functions Mj, for 
j = 2, . . . , n, such that 



F(w, w, t, v) — K t (w, w, t, 0, u) = 0. 



With such (13.81) simplifies to 



w = eM 1 (w, w, v) + e 2 H(w, w, t, e, v). 



(3.11) 




and so (13. lip becomes 



w = ewLi(ww) + eW(w, w, t, e, v) 



(3.12) 



where 



it 




(3.13) 



Differentiating the polar coordinates w = pe l ^ ^ yields 
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(3.14) 



But 

tie*®-*) = (ewL 1 (ww)+eW(w,w,t,e,v))eW-V 

= epL, (p 2 ) + ee^-VWipe'W-V , pe 1 ^ , t, e, v) 

and so 

p = eRl(p) +eR(t,i(>,p,e,u) 
i) = l + e^ (p) + e^(t,jp,p,e,u), 

where Rl(p) = pRe [Li(p 2 )], ^ (p) = -Im [£i(p 2 )] and 

R{t, i/j, p, e, v) = Re [^"'^(pe-^^, pe*^, t, e, i/)] 

1 (3 15) 

#(t,^,p,e, i/) = --Im ^-^^(pe-^-'^pe^-^^^, i/)] . 

Note that R, ^ G SPfTlSP^ and that \1/ is not defined at p = 0. We now give sufficient 
conditions for the existence of an integral manifold in ( 13.141) . 

Theorem 3.2 Assume that R and , as defined in {\3.15\i . are C l on intervals away 
from p = and that the averaged equation 

P = <(p) (3-16) 

has an equilibrium pi > with D p Rl(pi) =71 7^ 0. If the parameters are small enough 
to satisfy the conditions outlined in the proof below, then Qff.i^P has an invariant torus 
S e>v , whose stability is exactly determined by the sign of e-ji. 

Proof: By the implicit function theorem, there is a neighbourhood 

n 

U = (-e*,e*) x JJ(-i^>, fj>) 

in parameter space and a continuous function p : U — > M + such that p(0, 0) = pi, 
ei^(p(e,z/)) =0 and D p B^(p(e,u)) = % ;U ^ 0, 
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where 7^71 > for all (e, v) G U, i.e. the stability of the equilibria p(e, v) is the same 
as that of p for all (e, v) G U . 

When e = 0, the phase space of (13.141) is foliated by invariant tori and so, from 
now on, we will assume that e 7^ 0. Consider the change of variables p = p(e, v) + x in 
03.14]) . Differentiating the new coordinates, we get x = p and the equivalent system 

x = ei?o(p(e, v) + x) + eR(t, ip, p(e, v) + x, e, v) 

ijj = l + e* (p(e, u)+x) + eV(t, ip, p(e, 1/) + x, e, 1/). 

By Taylor's theorem, there are continuously different iable functions B\ and B 2 such 
that 

Rl(p(e, v) + x) = Rl(p(e, u)) + D x R*(p(e, u))x + B x {x, e, v)x 2 
* (p(e, 1/) + x) = * (p(e, 1/)) + B 2 (x, e, 

Since i?J(p(e, z^)) = and D x R\(p(e, v)) = -y e>u , we obtain the new system 

i: = e"f eiU x + eA(t, 0, x, e, z/) 
^ = d(e, v) + 0(t, 0, x, e, z/), 

where 

A(£, -0, x, e, 1/) = i?i(x, e, z/)x 2 + i?(£, ^, p(e, z/) + x, e, 1/) 

6(t, ^, x, e, v) = eB 2 (x, e, z/)x + e\P(t, ■0, p(e, ^) + a;, e, ^) 

d(e,i/) = l + eV (p(e,v)) 

are at least C 1 by hypothesis. 

Let /7 + = {<, EU : q > for alH = 1, . . . , n) and (e , z/ ) G U + . Define 

n 

So = [-e , e ] x JJ[-^bj,^bj]. 

5=2 

As pi > and p(e, 1/) is continuous on U, it is possible to chose (e , vq) in such a way 
that 

(T = min {p(e, z/)} - |p x > 0. 

(e,v)GSo 



(3.17) 
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If x > —cr , then p = p(e, v) + x > p(e, v) — a > |p x for all (e, v) £ So. In that case, 
B and A are continuously differentiable, as R and \l/ are continuously differentiable in 
p on [±px, oo). Note further that 6, A £ tyf n qjj* 

Set S = 1 x R x {0} x So, and S CT = R x R x [-a, a] x £„. Then 

1. 6 and A are bounded by a function S(e, v) = 0(e, v 2 , ■ ■ ■ , u n ) over S (see ap- 
pendix for details), and 

2. for all < a < cr , O and A are Lipschitz in Hale's sense (with Lipschitz con- 
stants 6(e, v, S) = 0(e, v 2 , . . . , v n , 5) and rj(e,u,5) = 0(e, v 2 , . . . , v n , S), respec- 
tively) over Eq. (see appendix for details). 

Accordingly, theorem 12.11 can be applied to show there is a neighbourhood Si C Sq of 
the origin in parameter space for which f)3.14p (since it is equivalent to (13.171) ) has an 
invariant torus 7~ e u when (e, v) £ Si. Furthermore, the stability of T € V is the same as 
that of the hyperbolic equilibrium p(e, u), which is given by 671. □ 

The invariant torus 7^„ appearing in the proof of Theorem 13.21 can be parameterized 
by a relation of the form x = T £)i/ (0i, 9 2 ), where 61,82 £ S 1 . Let 

(£,„) = ^ J^J\A^ d9 lC l9 2 (3.18) 

denote the center of T e u , and let £ e>v be the corresponding epicyclic manifold of (13.4)) . 
in which all solutions are epicycles when projected upon the z— plane. 
Define the average value 

+ e«((p(e, 1/) + <f^))e-^- i) , c.c, t, e, z/)) #tZt. (3.19) 

If 7^„ is stable (in the sense of theorem 12.11) . we shall say that [£ € ,v]d is the center of 
drifting of E^ v . 

Theorem 3.3 Suppose the hypotheses of theorem \3.^\ are satisfied. Then there exists 
a wedge-shaped region near X = of the form 

V>! = {(A 1; . . . , X n ) £ R n : |A 3 -| < VxjlXxl, V hj > 0, for j ^ 1 and X 1 near } 
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such that for all 7^ A G V\, <\2.2ty has an epicycle manifold £\, with [£l}r> near, but 
generically not at, £1. Furthermore, [£j[]d a center of drifting when A171 < 0. 

Proof: According to theorem 13.21 there are constants ei, 1/1,2, ■ • ■ , vi >n > and a 
neighbourhood 

n 

Si = [-ei, ei] x JJ[-z/ij, 

such that (I3.14p has an integral manifold £ eiV whenever (e, z/) G Si. 
For j 7^ 1, set Ai = e 7^ 0, Aj = uje and Vy = z/y. Then A G Vi as 

|Aj| < ■ [Aij < Vi^jAil for j 7^ 1, 

and f)2.2p has an integral manifold £\ = £i—ie lt v+S eiU . Furthermore, [^]d = £i+0(\i) 
and so [^]d 7^ f° r a generic 7^ A G Vi. The conclusion on the stability of £\ then 
follows directly from theorem 13.21 □ 

Remark 3.4 1. These isolated epicycle manifolds need not in general be unique for 
a given A G Vi as Rq{p) = may have any number of hyperbolic solutions. 

2. In generic semi-flows, all that can be said with certainty from the analysis when 
the parameter values stray outside of Vi is that the epicycle manifolds in (12.21) 
drift away from £1, which cannot then be a center of drifting. This is not unlike the 
situation with regards to spiral anchoring [7]. Richer dynamics and interactions 
with rotating waves can also take place; for instance in [5], we gave an example 
in which the epicyclic manifold collapses at a saddle-node bifurcation of rotating 
waves. 

3. Note that the actual parameter region in which epicyclic drifting is observed may 
be much larger than V\\ however, our local analysis cannot be used to obtain 
global results. 

The preceding results have been achieved by considering (12. 2p under a co-rotating frame 
of reference around £1. Of course, since the choice for £1 was arbitrary, corresponding 
results must also be achieved, in exactly the same manner, when the viewpoint shifts 
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to another Indeed, for j — 1, . . . , n, define the average functions 

M j {w,W) = — [ e*H j (we- a ,we it ,0,0)dt; 
2tt Jo 

as before, each M J is S 1 — equivariant and so there are continuous functions Lj : K — »• C 
such that Mi(w,w) = wLj(ww). We will call 



^ (p) = pRe[L J (p 2 
the epicycle functions of (12. 2p . 

Corollary 3.5 Lei £ {1, . . . , n}. If p* > suc/i £/iat 

i?o(P*) = and D p R k {p*) = 7* ^ 0, 

£/ien i/iere exists a wedge-shaped region near A = of the form 

V k = {(Ai, . . . , A n ) £ R n : |Aj-| < V kJ \\ k \, V kJ > 0, for j ^ k and A fc near } 

such that for all 7^ A £ V k , §£M has an epicycle manifold £\, with [£x]d near, but 
generically not at, £ k . Furthermore, [££]d is a center of drifting when A&7* < 0. 

Proof: The epicycle function is exactly the function that would appear in (13.161) 
had the preceding analysis been done around £ k . Theorems 13.21 and 13.31 can then be 
applied directly to obtain the desired result. □ 

Clearly, the remarks appearing after the proof of theorem 13.31 still hold. There is 
one last statement to be made concerning epicycle manifolds: theorem 13.51 only gives 
sufficient conditions for their existence in (12.21) . In section [5], we have provided an 
example that shows that they are not, in fact, necessary conditions. 

4 Epicyclic Drifting For Combined RSB-TSB Terms 

In this section, we investigate another way in which the Euclidean symmetry can be 
broken: by combining rotational and translational symmetry breaking. In effect, we 
are lifting the restriction (3 = 0, with n = 1 in ( 12.2ft . 
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It turns out that the value of f plays a crucial role in the analysis: the cases f = 1 
and j* > 1 are essentially different. The general lines are very similar to those of the 
preceding section, as such, the proofs are omitted in order to avoid tedious repetitions. 
In either case, however, me assume without loss of generality that £1 = 0. 

4.1 The Case f = 1 

Let F G : R x M 2 -> C be defined by 

F Q {m = e*\-iv + pY, ¥^TrT > 
L i[m + 1) J 

where the g m (P) are the Fourier coefficients of G 6 ^P^. Set 2 = p — F G . Then, (12.21) 
rewrites as 

i = Pg_ x {p) + Pe u H((z + F G (t, P)e~ u , c.c, X 1 ), (4.1) 

where c.c. represents throughout the complex conjugate of the preceding term. Gener- 
ically, g-i(0) ^ 0. Set e = /3, v — \\ and e = eAi. Then (14.11) transforms to 

z = L>H*(ze- tt ,ze a 7 t 7 e,L>), (4.2) 

where 

H*(w,w,t,e,v) = eg_ l {ev) + e it H{w + F G {t,ev)e~ i \c.c.,v) (4.3) 

is 2/T— periodic in t, smooth and uniformly bounded in w. Consider the near-identity 
change of variables 

z = w + vg(w,w,t,e,u) (4.4) 

where g G ^ n is continuous in all of its variables and to be determined later. This 
change of variables transform (14.21) into the equivalent system 

w = v (H*(we~ l \we l \t,e,ti) - g t (w,w,t,e,0)) + v 2 H*{w,w,t,e,v), (4.5) 
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where G ^ 27T is bounded and continuous in all its variables. Denote the average 
value of H*(we~ lt , we lt , t, e, 0) by 

i r 2w 

MJw,w,e) = — / HJwe~ lt ,we lt ,t,e,0)dt. 
27r Jo 



Since 



H.iwe-'^we*, t, e, 0) = eg^(0) + e lt H ({w + F G (t, 0))e" i ', c. c, 0) 

= %-i(0) + e u H(we~ lt - iv,we u + iv, 0), 

we have 

1 f 27T 

M*(w,w,i) = %-i(0) + — / e lt H(we~ lt - iv,we lt + iv, 0) dt. 

2tt Jo 

Then 

H*(we~~ u , ^,£,6,0) = M*{w,w,e) + F*(w,w,t,e), 
where G ^P^ 71 " is uniformly continuous and 

F(w,w,t,e)dt = 0. (4.6) 

o 

Let be an antiderivative of F* with respect to t. Then £> G ^ 77 by (14. 6 p and 

F*(w, 17, t, e) — W, t, e, 0) = 0. 
With such a (14.51) simplifies to 

w = vM*(w,w,e) + v 2 11.*(w,w,t,e,h>). (4.7) 

As M,,(u>, 17, 0) is also S 1 — equivariant, there is a continuous function L* : M. — > C such 
that M*(ty,iU, 0) = wL*(ww), and so (14.71) becomes 

li; = uwL^ww) + vW^w, w, t, i, v), (4.8) 

where 

W*(w,TD, t,e, v) = eg_i(0) + vH.*{w,w,t,e,v). (4.9) 
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Differentiating the polar coordinates w = pe l ^ *-* and substituting in 14.81 yields 

p = vR* (p) + uR(t, tp, p, e, v) ^ 
ip = l + W*(p) + rp, p, e, i/), 

where i?*(p) = pRe [L*(p 2 )], *»(p) = -Im[L,(p 2 )], i2, e qjf n ^ and * is not 
continuous at p = 0. We now provide sufficient conditions for the existence of an 
epicycle manifold in (12.21) . 

Theorem 4.1 Assume that R and , as defined in ( 4-10 ), are C 1 on intervals away 
from p = and that the averaged equation 

p = eR*{p) (4.11) 

has an equilibrium po > with D p R^(p ) = 7 ^ 0. If the parameters are small enough, 
there exists a wedge-shaped region near ((3, Ai) = (0, 0) of the form 

V = {{13, Ai) e M 2 : \(3\ < i^Ail, K > 0, Ai near } 

such that for all (/3, Ai) G V, /3 ^ 0, /ias an epicycle manifold Ql Xi , with [Ql a Jd 
near, frirf generically not at, the origin. Furthermore, [Q\ a Jd «s a center of drifting 
when Ai7o < 0. 

The remarks after theorem 13.31 still hold after having been suitably modified. 



4.2 The Case f > 1 

The case f > 1 is handled slightly differently. Let Cjj(C) and Cjj(C) be the spaces of 
continuous and continuously differentiable functions from M to C, respectively. Then 

Cl /f ={f:fe ^? lf H C°(C)} and = {/:/£ PI Cfc(C)} 

are Banach spaces when endowed with the respective norms 

|H|o = sup{|u(x)| : x G [0, 2n/f}} and | \u\ \i = \ \u\ | + | \u'\ | , 
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and the linear operator y : C\^,* — > C^if defined by y(u) = iu + v! is bounded, 
invertible and has bounded inverse. 



Define the nonlinear operator 7i G '■ ^Iw/j* x ^ 2 ~^ by 



n G (u, (3, Ai) = y(u) -\ lH [ u -tv + pYl T(mf + iy ac -' Al ) ' (4 ' 12) 



where the g m (P) are as they were in the previous section. But TC G (0, 0,0) = and 
Di?i G (0, 0, 0) = i ^ and so, by the implicit function theorem, there is a neighbour- 
hood M of the origin in M 2 and a unique smooth function U : K 2 — > C^A* satisfying 
U(0, 0) = and H G (U((3, Ai),/3, Ai) = for all (/?, A x ) G A/". 

Let F G : R x M 2 -> C be defined by 

r n ( FC\p lrn 3* t "i 

F G (t, (3, Ai) = e« [ - iv + /3 Zll T + Ax)(t) . (4.13) 

Then Ai))(t) = \\H (F G (t, (3, Ai)e _it , c.c, Ai) , and, upon setting z = p - F G , 

( 12.21) rewrites as 

i = X ie u [H((z + F G (t, (3, X^e'*, c.c, X 1 ) - H(F G (t, (3, Ai)e" lf , c.c, Ai))] , 
which reduces to 

z = X 1 e it H(ze- it ,ze it ,t,/3,X 1 ), (4.14) 

where 

H(w, w, t, (3, Ai) = H{w+F G (t, [3, X x )e~ i \ c.c, X x )-H{F G {t, f3, Ai)e" lt , c.c, Ai) (4.15) 
is 27r/j*— periodic in i. Then, ( 14.141) becomes 

p = \ 1 R*(p)+\ 1 R(t,il;,p,P,\i) 

(.4.1b) 

V> = 1 + Ai* *(p) + Ai^(t, ^, P, V, Ai), 



where R*(p) = p Re [L*(p 2 )], \l/*(p) = — Im [L*(p 2 )] for some continuous function 

r n <p 2 - 
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C, * G *p 27r n $8* and * is not continuous at p = 0. 



Theorem 4.2 Assume that R and , as defined in fl^.iffi ), are C 1 on intervals away 
from p = and that the averaged equation 

P = eR*{p) (4.17) 

has an equilibrium po > with D p R*(po) = 70 7^ 0. // the parameters are in a (small 
enough) deleted neighbourhood V J * of the origin, (\2.2l) has an epicyclic manifold Qi Ai , 
wzt/i a Jd = 0. Furthermore, the origin is a center of drifting when A170 < 0. 

Remark 4.3 1. The small term eg_i(0) in f ]4.3f) and the absence of a corresponding 
term in (14. 15f) are responsible for the different form of the regions V (in theorem 
14.111) and V 3 * (in theorem 14 .171) . as well as for the location of the center of drifting. 

2. The remarks made after theorem 13.31 still hold, when suitably modified. 

3. There are a lot of similarities between our analysis and the results obtained during 
the analysis of spiral anchoring in [6,7], such as the presence of wedge-shaped re- 
gions or the deleted neighbourhoods, depending on the nature of j*. In particular, 
one might hope that the epicyclic manifolds would possess Z.,*— spatio-temporal 
symmetry; however, this is not the the averaged system defined by (14.161) 
is generally only 2n/j*— periodic in t when R = and \l/ = 0. That being said, the 
epicycles themselves possess this symmetry in an appropriate co-rotating frame 
of reference. 

5 Epicyclic Drifting For General ESB Terms 

Lifting all restrictions on (3 and A in (12.21) . and combining the methods of the preceding 
section, we obtain the following general epicyclic drifting theorems for (12.21) . 

Theorem 5.1 Let n > 1 and k G {1, . . . ,n}. Given a hyperbolic equilibrium p* of an 
appropriate averaged equation p = (3Rq(p) with eigenvalue 7^ (derived as in section 3), 
there is a region in parameter space near (/3 = Ao, Ai, . . . , A n ) = of the form 

Vl = {(A , Ai, . . . , A n ) G R n+1 : |Aj| < V Jtk \Xk\, V jjk > 0, for j ^ k and A fc near } 
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when f = 1, or of the form 

V( = {((3, Ai, . . . , \ n ) € R n+1 : \P\ < /3b, |A,| < V j>k \X k \, V j>k > 0, for j ^ k and A fc near } 

for some constant (3q > 0, when f > 1, such that for all ^ A) G VjJ , wit/i the 
additional condition that (3^0 when f = 1, (\2.2t) has an epicyclic manifold Ep'* , with 
[Ep'^D near, but generically not at, Furthermore, [Epxln is a center of drifting 
when Xklt < 0- 

Since the hypotheses of this theorems are not generic, it is not clear that such integral 
epicyclic manifolds are common, and their existence must sometimes be inferred in the 
physical space, especially if they are repelling, such as appears to be the case in [25]. 



6 An Epicyclic Manifold in Physical Space 

In this section, we provide what we believe to be the first observed example of epicyclic 
drifting in a modified bidomain system. Our system is a TSB perturbation of the 
bidomain equations and parameter values found in [27]: 

- 3 \ CVF 



1 / U' J 

Ut = -A U ~Y 



Au 



yy 



2 + a + 



a 



1 + a(l - s) 
+ I 2 + a{\ - e) 



1 



a(l — e) 



HI) 



1 



1 



ail — e) 



v t = q(u + 5 — jv) + (f)(x — 35, y — 35) 



where q = 0.3, a = 1.0, e = 0.75, 5 = 0.8, 7 = 0.5 and 

(j)( Zh Z2 ) = -0.03 exp (-0.085(^ 2 + z$)) 



-1 



(6.1) 



00. 



The TSB term <j>(x — 35, y — 35) is uniformly bounded and goes to as || (x, ; 
Futhermore, it preserves rotations around the point (35,35). 

As our emphasis lies with qualitative observations rather than with precise nu- 
merical analysis, the numerical perspective is somewhat naive. The computations are 
carried out on a two-dimensional square domain [10, 60] 2 with 120 grid points to a side 
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I 1 1 1 1 1 1 

10 20 30 40 50 GO 70 

Figure 3: Tip path of the u component of two solutions of (16.11) . The solution in red is attracted 
by an epicyclic manifold. The solution in blue shows the effect of the boundary: the two solutions 
are clearly not of the same nature. Compare the epicyclic manifold with the image in figure [TJ 



and Neumann boundary condition, using a 5-point Laplacian, continuous linear finite 
elements on square meshes and the fully implicit second order Gear finite difference 
integrator. The tip path of the u component of two solutions are shown in figure 3. 

7 Summary and Concluding Remarks 

Recently, equivariant dynamical systems theory has been used to provide an approach 
to the study of spiral wave dynamics and bifurcations, in particular, it has provided 
mechanisms for such behaviour as spiral tip meandering and resonant growth [2,30-33, 
38], spiral anchoring/repelling and boundary drifting [19], which have been explained 
as consequences of forced Euclidean symmetry breaking. 

In this paper, we have used this model-independent approach to analyze epicyclic 
drifting of the spiral tip in media with several localized inhomogeneities, with or with- 
out anisotropy. The result of a simple numerical experiment is in agreement with our 
theoretical conclusions. The RSB terms are characterized by an integer j* > 1. It is 
important to note that, as of now, only the integers f = 2 (anisotropic cardiac tissue, 
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say) and f — 1 (an excitable medium that in which there is a directed current, such 
as a reaction-diffusion-advection system) have easy interpretations in the context of 
RSB. 

It should be recalled that our analysis rests on certain simplifying modeling assump- 
tions which may not be valid in some realistic physical systems: namely concerning 
the discrete nature of the inhomogeneities (i.e. finite number of inhomogeneity sites) 
and the hypothesis of local rotational symmetry of the individual inhomogeneities. In 
a realistic model of excitable media such as the bidomain model describing electri- 
cal conduction of cardiac tissue, with or without advection, actual inhomogeneities 
may lack this local rotational symmetry, or may even be distributed smoothly and 
non-symmetrically over the medium {i.e. an inhomogeneity field). 

Should the discrete localized inhomogeneities not possess circular symmetry, we be- 
lieve that our results would still describe the essential qualitative features of epicyclic 
drifting, even though our analysis does not technically apply, as epicyclic drifting is 
linked to hyperbolic fixed points. However, for a smoothly distributed non-symmetric 
inhomogeneity field, our techniques are unlikely to yield meaningful results. 

Finally, we would like to point out that our existence results do not provide a de- 
scription of the flow on the epicyclic manifolds: we gave two examples of center bundle 
systems with essentially different flows in [5]. In the first example, the flow is "ergodic" 
(see the first figure of this article): if the spiral tip lies in the epicyclic manifold, the 
tip path eventually fills the entire manifold. In the second example, the flow on the 
epicyclic manifold is dictated by the stability of rotating waves located on the manifold 
(see figure 4). 
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Figure 4: Projection (in the z— plane) of a stable epicycle manifold £ containing three rotating 
waves s, n and f. The arrows indicate the flow on the manifold. 



A Technical Results 

Proposition A.l The function M 1 (w,w,0) defined in (\3.9f} is S 1 — equivariant. 
Proof: Recall that H\ is S 1 — equivariant by construction. Then 



1 



2tt 



M 1 (we~ ie ,we ie ,0) = — / K (we~ ie e~ u , we i6 j\ t, 0, 0) dt 

C-27T 



= — [ * e it H 1 (we- it e- i9 , we u e 10 ,0,0) dt 
27r Jo 

= — [ " e u e ie H^we-", we?*, 0,0) dt 
27r Jo 

= e ie M 1 (w,w,0), 

that is, M 1 is S 1 — equivariant. □ 
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Proposition A. 2 Let all terms, variables and functions be as in Theorem \3.2i In 
particular, the functions R and \l/ are C 1 on intervals away from p = 0. Denote 



Note that these spaces, as well as the spaces E and E CT from Theorem \2.1[ are convex. 
The functions and A satisfy the following conditions: 

1. O and A are bounded by a function H(e, v) = 0(e, 1/2, ■ ■ ■ , v n ) over S , and 

2. for all < a < a , and A are Lipschitz in Hale's sense (with Lipschitz con- 
stants 9 [e,v, a) = 0(e, v-i, ■ ■ ■ , v n , a) and u{e,u,a) = 0(e, v-i, ■ ■ ■ , v n , a), respec- 
tively) overlu a . (p. [THjl 

Proof: Since 0, A G tyf* n ^P^ 71 ", we need only show that the first statement holds over Sq 
and the second over for all < a < gq. For j = 1, . . . , n, there are appropriate functions 
Rj, 9j G <p t 27r n q?^, C 1 on intervals away from p = 0, such that 



= [0, 2tt] x [0, 2tt] x {0} x S 
E; = [0, 2vr] x [0, 2tt] x [-a, a] x 5, 



n 



3=2 



(A.l) 



n 



R(t,ip,p,e,v) = eRi(t,i/j 



i=2 



according to (|3.13p . 



1. Over Eq, we have x = and so 



9(i, V, 0, e, 1/) = etf (i, ^, p(e, u), e, v) 
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(A.2) 



= e 2 ^i(t,i;,p(e,u),e,v) + e ^ VjVj(t, ip, p{e, v), e, v) 

J=2 



and 



hit, ip, 0, e, 1/) = ip, p{e, v), e, z/) 



n 



(A.3) 



= ei?i (t, ip,p(e, v),e,v) +^ (t,if), p(e,v), e, v) 

i=2 
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according to (jA.ip . 

For j = 1, . . . ,n, the continuous functions \Rj\ and le^jl reach their maximum Cj and 
Ej, respectively, on the compact set [0, 2ir] x {0} for j = 1, . . . , n. Then 

\e^j{t,t/j,p(e,u),e,u)\ < Ej and \Rj(t, ip, p(e, u), e, v)\ < Cj 

over Eq for j = 1, . . . , n. According to (|A.2[) and (|A.3p . 

n 

|0(i,^,O,e,i/)| < \e\E 1 + J2Wj\E j = Qi(e,v) 

7=2 
n 

|A(t,V,0, €,i/)| < |e|Ci+ J^I^ICj = Q 2 (e,z/) 

7=2 

over Eq. Set 

E(e,z/) = max{<2i(e, v),Q 2 {e,v)}. 

Then A and G are bounded by 3(e, v) = 0(e, z/ 2 , . . . , u n ) over Eq. 

2. Let (^1, x\), (ip2,%2) £Rx [—0", <r] for < a < <jq. By one of the mean value theorems, 
there exist points (ip* , x*), (^*, x*) G [0, 27r] x [— cr, er] on the line joining (ipi,x±) and 
(1^2, X2) such that 

|0(t,^i,a;i,e,i/) 
|A(t,^ 1; a;i,e,i/) 

where 

A(t,if),x,e,v) 

i=2 



- 0(i, ^2,^2, e, f)l = |0(t,^*,x*,e,i/)|[|^i - V2I + \xi -x 2 \] 

- A(t,if} 2 ,x 2 ,e,v)\ = |A(t,^*,a;*,e,z/)|[|^i -^ 2 | + |si - z 2 |], 

= D x Q(t, ip, x, e, v) + D^@(t, ip, x, e, z/) 

= xKq(x,€,u) + eKf(t,iJ},x,e,v) + UjKfjt, ip, x, e, v) 

7=2 

= D x A(t, V>, e, z/) + D^A(t, ip, x, e, z/) 

rt 

= xK^(x,e,u) +eKf t (t,ip,x,e,u) + ^ z^-iff (t,i/>, x,e,u), 
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where 

K$(x, e, v) = eD x B 2 (x, e, v) 

K^(x, e, v) = e (D X B 1 (x, e, u)x + 2B X (x, e, v)) 
Kf(t, ip, x, e, v) = B 2 (x, e,u) + e (D^^t, tp, p(e, v) + x, e, v) + A/,*i(i, tp, p(e, v) + x, e, v)) 
K^(t,ip,x,e,u) = e(D x R 1 (t,ip,p(e,u) +x,e,u) + D^ } R 1 (t, ip, p(e, v) + x,e,u)) 

and 

Kf(t,ip,x,e,v) = e(D x ^j(t,ip,p(e,u) + x,e,u) + D^j(t, ip, p(e, v) +x,e,v)) 
Kf-(t,tp,x,e,u) = e (D x Rj(t,ip, p(e, v) + x,e,u) + D^Rj(t,ip, p(e, v) + x,e,u)) , 

where tyj and Rj are as in (|A.1|) . Since G and A are continuously differentiable, and 
A are continuous on E£, as are Kj and for j = 0, . . . , n. 

In particular, the functions \K^\ and \Kj\ each reach their respective maximum kj 
and kj 1 on ££. for j = 0, . . . , n. Then, note that | 

n 

\@(t,r,x*,e,u)\ < \x*\\K^(x*,e,u)\ + \e\\Kf(t,r,x*,e,u)\+Y,W 3 \\Kf(t,r,x*,e,u)\ 

3=2 

n 

< \x*\k* + |e|*f + j>#? 

3=2 

n 

< ak* + |e|fcf + Wj\*>j = 

3=2 

and 

n 

|A(i,^*,a:*,e,z/)| < (x*,e,i/)| + |e||iiCf (t, V 1 *, x*, e, + ^ t/>*,ai*,e, 

3=2 

n 

< |a;*|ft£ + |e|fcf : + ^1^1*? 

3=2 

n 

< akg + |e|*f + Y Wj\k? = v{e, v, a). 

3=2 

In particular 

\&(t,ipi,xi,e,v) - Q(t,i() 2 ,x 2 ,e,v)\ < 0(e,u 2 , v n ,a)[\ipi - ip 2 \ + \x x - x 2 \] 
\A(t,ipi,xi,e,i/) - A(t,i() 2 ,x 2 ,e,v)\ < v(e,v 2 , . . . ,v ni a)[\il>i - ip 2 \ + \xi - x 2 \], 
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where 9{e, v, a), v(e, v, a) = 0(e, U2, ■ ■ ■ , v n , a). Hence and A are Lipschitz in Hale's 
sense. 

This completes the proof. □ 
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(a) fx = 0.001 (b) j.i - C 



